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Abstract 

^N . A measurement of the Lamb shift of 49,881.88(76) GHz in muonic hydrogen in conjunction with 

^ , theoretical estimates of the proton structure effects was recently used to deduce an accurate but 

bX)' rather small radius of the proton. Such an important shift in the understanding of fundamental 

•^ ■ values needs reconfirmation. Using a different approach with electromagnetic form factors of the 

^^ proton, we obtain a new expression for the transition energy, A = E2P — -Eg^ , in muonic 

,^ ' hydrogen and deduce a proton radius, Vp = 0.831 fm. 
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I. INTRODUCTION 

In a Quantum Field Theoretic (QFT) description of a muon-proton amplitude, the 
coupling of the photon to a point-like particle (vertex), ej^, has to be replaced by 
e(Fi(g^)7^ + F2{q^)a^yq^ /2mp) where q^ is the four momentum carried by the photon, 
(f = Q'ai^^ s-'^d Fi(q^) are the electromagnetic form- factors jll, 121] which encode the infor- 
mation on the structure of the proton. The inclusion of form factors introduces corrections 
to the binding energy in Coulomb bound systems such as the hydrogen atom. In fact a sys- 
tematic inclusion of the proton structure effects introduces corrections to other terms such 
as the fine and hyperfine splittings in the hydrogen atom 3j. Given the recent high precision 
of experimental measurements, the corrections due to the proton size become relevant. One 
such precision measurement using pulsed laser spectroscopy was recently performed 4] to 
measure the Lamb shift in muonic hydrogen (fi^p). Using this accurate measurement of the 
energy difference (Lamb shift) between the 25*1/2 and 2P3/2 states, a precise value for the 
transition energy, A = E2P — E^g = 206.2949(32) meV, in muonic hydrogen was given. 
Comparing the measured value with a model calculation of this difference which included 
the proton structure effects, a new value of the proton radius, namely, Vp = 0.84184(67) fm 
was published. It was further claimed that the new value implies that either the Rydberg 
constant has to be shifted by -110 kHz/c or the calculations of the QED effects in atomic 
hydrogen or muonic atom are insufficient. The uncertainties due to form factors have been 
discussed in [4]. Another uncertainty however seems to arise from the way in which one 
incorporates the finite size corrections (FSC) due to the structure of the proton. The FSC 

^ were based on methods given in p, |6|. Here we implement the FSC based on the Breit 
7|] and thus obtain a different form for the expression of the energy difference, 
namely, A(= E^p^^^ - ^Ij^/J = 209.16073 + 0.1139rp - 4.3029r2 + 0.02059r3 meV, as com- 
pared to, A^"*"'^^ = 209.9779(49) - 5.2262r2 + 0.0347r3 meV, obtained in ^]. Based on this 



m 
equation 



new expression and following the same procedure as in 



B, 



we extract a new radius of the 



proton. This value of Vp = 0.83112 fm based on the FSC from the Breit equation is not far 



B 



from that predicted in [^, namely rp = 0.84184(67). However, given the precision of the 
numbers under consideration, it is worth noting that the method used for proton structure 
corrections does introduce an uncertainty. 

Before proceeding to the formalism used in the present work, we note certain points which 



motivated the present work. The proton radius is one of the fundamental numbers in nuclear 
andparticle physics and is obviously relevant for atomic physics too. The radius estimated 
in \^ is smaller than the value obtained by some other methods (see [4] for references and 
discussions on this issue). This could mean that some corrections to the energies in the 
hydrogen atom are missing. If a shift in the understanding of the proton radius occurs, it 
is necessary to reconfirm the result with an independent calculation. This is particularly 
so in view of the model dependence entering the extraction of the radius. The nuclear 
physics inputs (here, the structure of the proton) are never as precise as their electroweak 
counterparts. In the present work, we try to provide this reconfirmation using a different 
theoretical approach for the proton structure corrections. Though we do agree qualitatively 
on a smaller radius, we notice that the accuracy gets blurred by the approach used for finite 
size corrections. 

II. BREIT POTENTIAL WITH FORM FACTORS 

All r dependent potentials in Quantum Field Theory (QFT) are obtained by Fourier 
transforming an elastic scattering amplitude suitably expanded in 1/c^ (see \^ for several 
examples). The Breit equation 8|, |9(| follows the very same principle for elastic e^fi^, e^e~ 
(positronium), e~p (hydrogen) and /i~p (muonic hydrogen) amplitudes. The one-photon 
exchange amplitude between the proton and the muon leads then to the Coulomb potential 
)lus the fine and hyperfine structure (hfs), the Darwin term and the retarded potentials 
8|, |9|. Here we use a modified Breit potential |3| for the /i~p system which includes the 
electromagnetic form factors of the proton. We present a calculation based on this potential 
with form factors to evaluate the transition energy. A, and hence the proton radius, Tp, 



as m 



a. 



The Fourier transform of the momentum space Breit potential in [3j] gives V = 
—a/r+6V{pf^,Pp, r), where SV{pf^,Pp, r) contains Coulomb, Darwin, fine and hyperfine terms 
with form factors. The expectation values of these terms give the energy corrections to the 
various terms. 

The standard Breit potential is normally written down either for point-like particles 
(with standard point-like vertices) or at zero momentum transfer at the vertices (i.e. -Fi(O) 
and ^2(0)), at least for the hadronic vertex. Therefore, extending it to include the full g^ 
dependence of the proton form factors (which encodes the finite size corrections) is just a 



straightforward and mild extension of the standard procedure. Thus one can obtain the 
Breit potential in momentum space for the e~p or fi~p system with the proton structure 
effects fully included and is given as (Eq. (21) in [Sj): 
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where, X = e or fi. The Fourier transform of this potential results in a space dependent 
potential suitable for calculations of the corrections to the atomic energy levels via time- 
independent perturbation theory. 

The two form factors, Fi{q^) and F2{q^) of the proton are connected to the Sachs form 
factors, GE{<f) and GM{(f')- The latter can be interpreted in the proton rest frame to be 
Fourier transforms of the charge and magnetization distributions in the proton. They can 
be approximated fairly well by a dipole form [l| (which is suitable for analytic calculations) 
as, G'D(q') = 1/(1 + qV^')' ~ G^q^) ^ G^,(q2)/^p, where (1 + Kp) = fip = 2.793 is 
the proton's magnetic moment and m the dipole parameter. We use the standard non- 
relativistic approximation g^ ^ — q^ to obtain the dipole form of Ff (q^) and -^^(q^) (see 
Eq.(30) in J3|). The parameter m is eventually related to the proton radius and taken as a 
free parameter to fix the radius of the proton (to be discussed below). 



III. PROTON STRUCTURE CORRECTIONS TO ENERGIES 



We follow a procedure similar to that in [4] to obtain the proton radius with the difference 

that the proton structure corrections are included via the Breit potential method. The 

f] 
Breit potential V{r) as explained above consists of a series of terms [8| with the leading 



term being the one due to the Coulomb potential. Thus, the expectation values of V{r) = 
Vc + Vd + Vhfs + Vso + ■■■ give rise to the various energy terms including the standard 
Coulomb, Darwin, fine structure, hyperfine structure etc. Below, we shall give the expression 
for some of the energies (using the Fourier transform of the momentum space potential with 
form factors [3]) which are relevant for the calculation of the radius. For details of the full 
potential in coordinate space, we refer the reader to |3|. 

A. Coulomb term 

The first term, namely, the Coulomb potential, V^f ^, with form factors is given as Vq^ = 



Vc + AVc, i.e 
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where k = 2mp/m with mp being the mass of the proton and m the parameter entering the 
dipole form factor. Note that the above expression has been obtained under the condition 
that 2inp ^ m and hence does not generate any singularities (which as such should not be 
expected too since the form factors are not singular). A small explanation regarding this 
and the other analytic expressions which follow is in order here. The potential in coordinate 
space is obtained by Fourier transforming the potential in momentum space. If we consider 
the Coulomb term (one with F^ Fl{Q'^)/Q'^ in Eq.([l]), where we write Q^ = q^), Fourier 
transform it and replace Fl{Q'^) using the dipole form factor, we obtain 
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We now use a partial fraction expansion to write the terms in the integrands. As an example, 
consider the second term in the square bracket of Eq. ([3]) . We can see that 
111 



(m^ + (^2)2(4^2 _^ g2) (4m2 - m2)2(4m2 + (^2) (Arn?p - rn?y {w? + Q"^ 
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(4m2 — m?){m? + Q'^Y 
assuming of course that Am? ^ rn^ . The integrals in (|3]) can be performed analytically and 
lead to Eq. ([2]) which contains terms of the type 

1 1 1 
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with k = 2mp/m. Should it happen that Arrip = m?, then l/[{w? + Q^)^(4mp + Q"^)] = 
[l/(m^ + Q'^Y] and the integral J^ dQ [Q sm{Qr) / {m'^ + Q"^)^] can be evaluated by differ- 
entiating twice the integral Ji = /q°° dQ [Q sin((5r)/(m^ + Q"^)] = {ji /2) e'"^^ with respect to 
m. However, the case of Arrip = im? seems unrealistic and is not considered in the present 
work. 

The correction to the Coulomb energy due to form factors is then found as, 
A£^f7 = {nlin\^Vc\nlm) , which with AV^ being only a function of r becomes, ^Ec = 
Jl^ AVc\Rni{r)\^r'^dr (with Rniir) being the unperturbed hydrogen atom radial function as 
found in books). The correction to the Coulomb energy for any n, / is thus given by a lengthy 
expression involving hypergeometric functions (as explained in the appendix). Performing a 
series expansion of the hypergeometric function and truncating it at large orders of the fine 
structure constant a, we evaluate the correction to the energy in terms of the dipole param- 
eter m. In order to rewrite the corrections in terms of the proton radius r^, we perform an 
expansion of the type G^(g^) = 1 — 2(^ jw? + ... of the Sachs dipole form factor and compare 
it with the standard G\{cj^^ = 1— < r^ > g^/G + .... One can then write < r^ >= \1jvp? 
and convert the corrections for energies given in terms of vn^ to those in terms of < r^ >. 
For convenience, in what follows, we shall denote < r^ >i/2= 7-^^ < 7-2 >= j-2^ < 7.2 >3/2_ ^z 
etc. Thus the corrections to the Coulomb terms of the 251/2 and 2Pi/2 levels used in the 
present work can be rewritten in terms of the proton radius r^ as follows: 



^251/2 m^CL' , ^ ^ , \ 2 "m^a^ , ^ . „^2 



A^r ^^(A,-A, + A,) rl + ;^(-2A, + 2^^ _ ^A^^ (4) 



with, Ai = 1 + [V(1 - k^fl ^2 = ^^p/[k^{l - k^)\ ^3 = 1 + [V(1 - k^)] and k = 2mjm. 

where, A = 3 + {nj{l - k^f] + [2^^/(1 - P)] - [l/{k\l - P)^)] and 5 = 7 + [2^^/(1 - 
k'^Y] + [5'tp/(l — k"^)] — [2/(A;^(l — k'^Y), makes a negligible contribution to the finite size 



a 



corrections. The Coulomb terms of the present work and those used in |^ will be compared 
in detail in section IV. 



B. Darwin term 



Let us first look at the Darwin term in the standard Breit equation. The potential in 
coordinate space is given as: 
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With the inclusion of the electromagnetic form factors, 
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which in the static limit, i.e. taking Fi2{q = 0) reduces to 
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Here X = e or fi depending on if we are considering the usual hydrogen atom or muonic 
hydrogen atom. In the above. 
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Taking the expectation values of these potentials using first order time-independent per- 
turbation theory, the corrections to the energies corresponding to the Darwin term can be 
found. For example, the Darwin term without form factors and for / = is given as, 
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and one with the inclusion of q^ dependent proton form factors is 
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where Goi{n,l) and G 02(1^, I) are lengthy expressions involving the hypergeometric func- 
tions, 2-^1- The Bohr radius, a^ = l/{m,ra), where m,r is the reduced mass. If we restrict to 
form factors taken at g^ = 0, then. 
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In order to get a better insight into the expressions, we rather replace the hypergeometric 
function 2-F1 by its series expansion 2-^1(0, b;c; z) = 1 + {ab/c){z/l\) + [a(a+l)6(6+l)]/[c(c+ 

1)] (2:^/2!) + in the expressions for energies and truncate the series at large orders of 

the fine structure constant a. Substituting for the dipole parameter by m^ = 12/rp, the 
corrections to the energies corresponding to the Darwin terms with n = 2 and / = 0, 1 are 
given as 



E 



D 



{2S 



1/2J 



and 



E 



D 



(2P, 



l/2j 



S 4 



+ 



S 4 

^2^ 



'1 + 2kx)+ 2 
■mi 



" 1 + 2kM + f (1 + 2kx)A^ + ^5i 1 rp(12) 



mt 



2kx)A2 



m 



X 



mt 



Bo\ri 



(1 + 2kx)A^ + 



ni 



X 



mt 



BArt + ... 



771 



1 + 2kx)D2 + ^E2 r;+ (1 + 2kx)D3 



771 



X 



771, 



7nt 



E,\rl + ... 



(13) 



The coefficients, Ai, Bi etc in the above equations depend on Kp and k = 2mp/7n. The 
correction for the 2P state turns out to be negligibly small and is not of much relevance for 
the present work. The coefficients which contribute to the 25* energy correction up to order 



Tp are given as 

A, 

A2 

Bi 

B2 
Here, a^ = l/(mj,Q;). 
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C. Fine structure 



As mentioned in the beginning, the expectation values of the various terms in the potential 
in q-space, Eq. ([1]), give rise to the Coulomb, Darwin, retarded potential, fine and hyperfine 
structures in the hydrogen atom. The fifth and the ninth spin-orbit terms in this equation 



are the ones corresponding to the fine structure. A Fourier transform of this potential leads 
to the following potential for the fine structure with form factors: 
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Using first order time-independent perturbation theory for evaluating the expectation value 
of the above potential, the fine structure energy term with the effect of form factors is given 

as, 
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where, Gpsi^, I) is given by a lengthy expression which contains the Gamma functions and 
hypergeometric functions. The electron (or muon) total angular momentum j = I + Sx, 
where I is the orbital angular momentum and sx the spin (with X = e or fi). The fine 



structure levels relevant for the present work are the 2P3/2 and 2Pi/2 levels which appear 

,2P3/, _2P, 

'FS 



in the difference, AE'i^q^'^ = Epg^^ — Ep^^^. The proton structure corrections to the 2P 



levels are always much smaller than those to the 25" levels, however, for completeness we 
give these expressions below. 
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with, Qr = l/(r7ir«) ^s before. The energy of the 2P3/2 level with form factor corrections is 
given by the relation Epg^'^ = —{l/2)Epg^^ . 



D. Hyperfine structure 



The hyperfine potential and the corresponding correction to the energy due to form 
factors in the case of muonic hydrogen has been discussed in detail in Ref. [3|. Here we give 



a brief description of this calculation for completeness. The 25* and 2P hyperfine energies 
relevant for this work are evaluated using the hyperfine potential, 
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The calculation of energies for states with / = and / 7^ is done separately 
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in section IV that the hyperfine levels of relevance are E^fY^ , Ef^^J^ , Ef^/J^ and Ef^ ^^^ 
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The complete expressions for any n and / can be found in [a]. For the 25*1/2 case, 
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""^f' ~ m,mp 12 y v^a/^+12af^ 123/2af^ + -yl ' ^'^^ 

with Or = l/{mra). Ej^j-]^^ = —SEf^j^l''^ and (l/4)A£'^jg = Ej^j]!'^ . Note the presence of 
the term linear in r^. The sum of this correction and one arising from the Darwin term will 
contribute a small term linear in r^ to the final expression for the transition energy, A. Such 
a term linear in Vp does not exist in the calculation of Pohl et al. [4]. The energies of the P 
levels are given as. 
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with Ci, C2, C3 and C4 the same as in the '2Py2 case. 

The contribution of the P levels to the finite size corrections are once again small. 



IV. RADIUS OF THE PROTON 

In the next step, we evaluate the radius by calculating the energy corrections to expres- 
sions dependent on the proton radius r^. Following a similar procedure as in [j] conceptually, 
we perform a fit using the experimental value of the transition energy, A = E2P, — E^g 
and the theoretical expressions for the corresponding energies (with finite size corrections) 
to determine the radiiis, \ 
we need to find as in J^], 
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to determine the radiiis, r^. In order to evaluate the energy difference, A = E2p^ — -^25, 
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Here / = j + Sp is the total angular momentum of the muon proton system with the muon 
total angular momentum, j = I + s^. The fine and hyperfine splittings respectively are. 
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A_E^j^ is the hyperfine splitting of the 25" levels, i.e., ^.E^j^ = Ef^^]!'^ — -E'/^j^^ . Such a fit 
in our case leads to the proton radius Tp = 0.83112 fm. The correction to the Coulomb term 
evaluated from Eq. (jl]) gives IS.Eq^ = 4.30248rp — 0.020585rpmeV, which forms the major 
contribution to the Vp and r^ terms to be discussed below. 



A. Coefficients of the r^ and rJ) terms 



The value -5.2262 r±+ 0.0347 rf, meV used in 
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equation taken from 10|] which is based on Friar's formalism [6|: 
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In this formalism, the correction to the 1/r Coulomb potential takes the form 

1 1 
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where p{s) is the charge distribution in the proton. Using perturbation theory, the correction 
to the energy level is evaluated in [6|] as (see Eq. (43a) in [6|]), 
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In contrast to our formalism where we consider the full r dependent wave function and 
perform the integration, Ref. |6| considers the wave function only at r = 0, as a result of 
which one is left with the integral < r^ > = / d^r p{r) r^ in the energy correction. First we 
discuss the < r^ > dependent term and then go over to the term with < r^ >(2)- 

The energy corrections and hence the coefficients of the r^ (writing < r^ >= r^ as in j^) 
term as calculated from Friar's formalism will clearly differ from the present work due to 
the difference in Eqs fl2^ and ([2]) for the two potentials. The finite size correction to the 
to Coulomb energy of the present work has been evaluated using the potential which is a 
Fourier transform of the first term in Eq. ([T]) , namely. 
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If instead of taking just the Coulomb potential in momentum space, we decide to club one 
of the Darwin terms with it, we find 



jjnewdef / \ 



1 



(26) 



q- 



A Fourier transform of this potential with a dipole form factor leads to the following potential 
in coordinate space: 



•trnewdef 



a 
r 



1-e" 



1 + 






(27) 



The energy correction evaluated for the 25* state (i.e. n=2, /=0) using AVc = V^^^ ^ 
i-a/r) in AEc = j;^ AVc\Rni{r)\^r^dr is then 

2 



mar 



mar 



1 + 3F(-2, 2; 3; 2/(1 + ma,))) + 



1 + m.ar 



1 + m,a. 

l + 3F(-2, 3; 3; 2/(1 + ma,))) 



(28) 
(29) 
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Expressing the hypergeometric function as a series expansion, truncating it at high orders 
in a and substituting ri = 12/iri^ leads to 

/.j^newdef _ Ci m^ In OmTV 3 \ /„„x 



The first term in Eq.( l30|) is exactly equal to that in (!23|) for the 25" state. This means 
that in principle, taking one of the Darwin terms together with the Coulomb one in the 
Breit potential leads to an exact agreement with the coefficient of the r^ term used in 4|. 
However, it also means that there are cancellations due to the other Darwin terms which 
eventually leaves us with a smaller coefficient of the r^ term than that used in ^. 

The second term in Eq.f l23l) is the third moment of the convoluted proton charge density 
and is defined as, 

< r^ >(2) = J d^r r^ p(2) (r) (31) 

where the convoluted charge density is given by 



P(2) 



J dhp,h{\z-r\)p,h{z). (32) 



Inserting the Fourier transform of the Sachs electric form factor Ge{Q'^) for p^h-, one finds 
that 

This third moment comes about due to the use of the smeared Coulomb potential (Eq. ( I24l) ) 
along with the perturbative expansion of the hydrogen wave functions. Using a dipole form 
factor Ge, one can show that < r^ >!?)= 35v^ < P >^/'^. This substitution in Eq. ( 123|) 
indeed leads to the factor 0.0347 r^ in |^. Such a correction of order a^ which is proportional 
to the proton form factor squared is not included in the present work. Including a term with 
the proton form factor squared in the Breit potential would correspond to a two photon 
exchange diagram which would obviously lead to a correction one order higher in a and 
hence would be quite small. 

B. Darwin contribution to the Lamb shift and the fine and hyperfine corrections 

The main contribution to the Lamb shift AE^s comes from QED, namely, the Uehling 
potential. In the Nature paper [^ using other methods |5|, |6| for including the proton 
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structure corrections, the authors found 



^^Nature ^ 206.0573(45) - 5.2262r2 + 0.0347rJ meV, (34) 

where the first term includes the relativistic one loop vacuum polarization for a point nucleus 
plus other QED corrections and the second and third terms are due to finite size corrections 
(FSC). To evaluate AE^s with FSC with the Breit potential approach, apart from AEq^ = 
4.30248rp — 0.020585rpmeV, we need the corrections AEq^Ii and the 25" and IP corrections 
to the Darwin terms too. The structure corrections to the Darwin and Coulomb terms of 
the IP levels are negligibly small. If we use Eqs ([9]) and ( ITTj) for the 251/2 level, we get, 
^noFi,2(g2) ^ 13.758591 me V and eJ'^'^^'=°^ = 14.4185121 meV. Starting with Eq.([IO]) and 
following the same procedure of expanding the hypergeometric functions etc, we find for the 
251/2 level, 

^Fi,2(g2) ^ 14418512 _ 0.0793 r^ + 0.0002613 r^ - 6.6 x 10"^ r J + . . . meV. (35) 

Note that the first term in the above equation is nothing but the Darwin term with form 
factors taken at c^ = 0. Since the calculation of the Lamb shift as used in [4] has been taken 
from a relativistic calculation, in principle we do not need to take into account the Darwin 
term. However, that calculation was done without including proton form factors and hence 
we must include the "correction" to the Darwin term due to form factors. Thus, we subtract 
^noFi,2(g2) ^ 13.758591 meV from Eq. ([35]) and use 

^E^Dallin = 0.64992 - 0.0793rp + 0.0002613r2 - 6.6 x 10"^ rj meV, (36) 

which is the correction to the Darwin term due to the proton form factors. Besides these, 

op op, 

we also have, AE^^^f = 3 x lO^^r^ - 1.37 x lO^V^ meV, AE^;/^„ = 9.08 x lO-^r^ - 4.4 x 



10 ^'^r'^ meV which are however extreme^ 



y small corrections. 



S with the value 206.0573 which is the sum of 24 

■ ^ n 



p 
To evaluate AEls, we start as in 

terms from the table given in the supplementary material of [4]. To this we add the FSC 

mentioned above for 2Si/2 and 2Pi/2 levels and find AEls = -^2Pi/2 ~ -^251/2 to be 



-') 



AEls = 205.40738 + 0.0793rp - 4.30274rJ + 0.020585r^ meV, (37) 

which should be compared with AE^g^'^^'^ mentioned in (1M|) above. Our result for AEls 
with FSC includes an additional term in r^ which is not present in AE^g^^^'^. The coefficients 
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of the Tp and r^ terms are however quite similar. The corrections to the P levels are really 
tiny as expected. 

Finally, the proton size corrections to the hyperfine and fine structure energy levels are 
evaluated as explained in the previous section. The fine and hyperfine structure terms in 
( 121 p and (12 2 p are expressed in terms of the proton radius as follows: 

AeJ5^' = 8.34678 - 4.26 x 10" VJ + 1.36 x lO^VJmeV (38) 



'^hfs — ^-^^-LZ, X. I o I /s xu I p -r o.-±u /S XU I p. 



AEt/J^ = 3.3912 - 1.787 x 10"V^ +5.45 x 10"Vi!meV, 



{l/4:)AEl% = 5.708 - 0.0347rp + O.OOOlrJ - 3.27 x lO^VJmeV. 

One can see that the FSC for fine structure are very small. The first term 8.34678 meV is 
exactly the same as the sum of first two terms given in Table I in the second reference of 



ll|. 



a 



In the evaluation of the proton radius in Pohl et al. [4j, the values of the hyperfine 
splittings were taken from (llj, where the FSC for the 25* level were evaluated using the 
Zemach method and those for the 2P case were not taken into account. Their FSC (taken 
from Table II of the first reference in [11|]) of order a^ and a^ sum to -0.1535 meV. This 
correction is obtained using a Zemach radius of Rz = 1.022 fm and leads to AE^fr^ = 22.8148 
meV. We note that using such a value of AE]ff^ means that the information about the radius 
of the proton (through the form factors or Zemach radius Rz = 1.022 fm) has already been 
included in the hyperfine energy which is later used to extract the radius of the proton. 
However, if we wish to remain within the spirit of the formalism used for AE^s, then we 
cannot take an input energy which already assumes a certain radius of the proton. Hence, 
instead of using just a number for AE^^^ (as done in j^) which already includes the FSC, 
we use the expressions given above in terms of Vp. 

Note that the above equations include the proton structure corrections only. Adding 
QED and other corrections to these levels as in Refs. [14,15] in [4] and following a procedure 



a 



similar to that in j^ (thus differing only in the inclusion of the FSC), we obtain 

A(= E(p,^^,^ - ^2^5^/2) = 209.16073 + 0.11388rp - 4.3029r2 + 0.020585r^ meV, (39) 

as compared to, 

^Nature ^ 209.9779(49) - 5.2262rJ + 0.0347r^ meV, (40) 
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obtained in J4|. 

Eq. (|39|) contains a term linear in Vp which is not present in the expression f l40|) of Pohl et 
al [^ . This term arises partly due to the proton structure correction to the Darwin term and 
partly due to the 25* hyperfine energy level. For the 25 hyperfine level, this is a correction 
of order a^ as can be seen from Eq. (1181) . In ^, such an order a^ correction of -0.1518 meV 
(see Table II in the first reference in ll| which was used in 4|) has been included directly 
as a number while taking into account the energy of the 25* hyperfine level (see also the 
discussion in the paragraph below Eq. (!38|) ). The remaining part of the Vp dependent term in 
our formalism arises from the form factor correction to the Darwin term which contributes 
to the expression of the Lamb shift in (1371) . In [41] the main contribution of 205.0282 meV 
to the Lamb shift arises from the relativistic one loop vacuum polarization which has been 
evaluated in [lOl] for a point nucleus. Hence, though the relativistic effects (represented by 
the Darwin term in our formalism) are taken into account, the finite size corrections to this 
term are missing in jj|. 

Taking the central value of the measured A = 206.2949(32) meV and replacing in the 
left hand side of Eq. (l39l) leads to Vp = 0.83112 fm. If we compare this radius with Vp = 
0.84184(67) fm (obtained from pUj) ) in [4] we see that an additional uncertainty arises due 
to the difference in the approaches used for including proton finite size effects. A detailed 
comparison of the present approach with those used in 4y can be found in section IV of Ref. 



To conclude, we can say that it is gratifying to see that in spite of the differences in the 
approaches for calculating the FSC, we obtain rp = 0.83112 fm which is not too far from the 



value of rp = 0.8418' 
different approaches 



(67) fm found in [^]. It is however important to note that there exist 



a, ISl ISl for evaluating the finite size effects in literature which can lead 
to different results (and an additional uncertainty in the extracted radius). One should be 
cautious not to overestimate the accuracy in calculating the radius of the proton. 



Appendix: Evaluation of energies 

The energies corresponding to the Coulomb, Darwin, fine structure and hyperfine struc- 
ture terms in the Breit equation are evaluated by taking the expectation value of the cor- 
responding potential using first order time-independent perturbation theory. In general, for 
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an operator A, the expectation value is, 

(A) = I r'^dr dO dd ^!,™, fr, 9, (/>) A^„,„. (r, 9 
where, 



" nlrrii V i^ i ' 



iv ^" ^-r /I 



vl/(r,^,0) = i?„Kr)l^ 



mi I 



with 



Rr,Ar) 



2 \3 (n-l-iy. 
nay 2n{n + 1)\ 



H 1/2 2r \ ' 



(A.l) 



(A.2) 



(A3) 



being the radial functions. For operators which do not depend on angles, the expectation 
value reduces to calculating 

1 «^i/2fn-/-7-l)\ (2j)! 



(A) 



/Ox 21+3 

\naj 



"v^' /2(n -l-j-l)\ 

2n22(«-'-i) j^^ \ n-l-j-1 J j\T{2l + j + 2) 



I" rfr Ae-2^/"V2'+2Ly +^) (4r/na) 



(A.4) 



The above integral can be evaluated analytically for the potentials considered in the present 
work. This amounts to e'^ 
^-mr 1^ etc. For example. 



work. This amounts to evaluating the expectation values of terms of the type e "^^, re '""\ 



IV 



1 



2/+3 

na) 2n22(«-'-i) 
r(2/ + 3) / na 



r(4/ + 3) V2 + mna 



h V n-l-j-1 )j\Ti2l+j + 2) 

21+3 / 4 \ 

F -2j, 21 + 3; 4/ + 3; — (A.5) 

V 2 + mna/ 



2 \2l+3 



na/ 2n22(«-'-i) 

r(2/ + 2) 



"d^i f2{n -l-j-l)\ r(4/ + 2j + 3) 



r(4/ + 3) \2 + mna 



,=0 V n-l-j-1 Jj\T{2l+j+2) 
'"'' )''^V(-2j,2/ + 2;4/ + 3;- ^ 



+ mna 



(A.6) 



1 



2 \2l+3 

na) 2n22(«-'-i) 

r(2/ + i; 



"d^i f2{n -l-j-l)\ r(4/ + 2j + 3) 



r(4/ + 3) \2 + mna 



,=0 V n-l-j-1 Jj\T{2l+j+2) 

'"'' Y^'f(-2j,21 + 1-41 + 3-—^ 
J \ 2^ mna, 



(A.7) 



2 \ 21+3 



1 



na/ 2n22(«-'-i) 



na 



r(2/) 

r(4/ + 3) V2 + mna 



E 

F f -2j, 2/; 4/ + 3 



2(n-/-j-l)\ r(4/ + 2j + 3) 
n-l-3-\ jj!r(2/+j+2) 

4 



2 + mna 



(A. 
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with the expressions being dependent on the Gamma functions and hypergeometric func- 
tions. The latter can be expressed as, 



V.c 2!c(c + 1) 



E 



(a)„(6)„2;" 



n=o lcj„ n\ 



(A.9) 



The expectation value of the Coulomb potential with form factors in Eq.([2]) for example 
leads to the Coulomb potential plus the finite size correction to the Coulomb energy for any 
ra, / which is given by 

/ 2 \ 2i+3 



AEcouiinJ) = a 



nttr 



1 



X 



"^Y2(n-/-j-l)\ r(4/ + 2j + 3) 

2n22("-'-i)r(4/ + 3) Pq \ n-l-j-1 ) j\T{2l+j + 2) 

21+2 / 4 



(l + 7r^)rP' + 2)(^^^) f(-2.2,.2;4,-.3;^^^) 



T3W^P'^^' 



na-r 



2 + rakna, 



2 + mnar , 

21+2 

F[ -2j,2/ + 2;4/ + 3; 



nar 



+— 1 + — ^ r 2 +3 , 

2 V l-Py ^ \2 + mnar 



21+3 



F( -2j,2Z + 3;4/ + 3 



2 + mknaj. 
4 



(A.IO) 



2 + mnar 

The energy corresponding to the Darwin term for example can be calculated by taking 
the expectation value of the potential in ([7]) and leads to Eq. fllU]) . namely, 



E 



D ' 



n,l) 



a 



2r7i^c2 n^af 



m 



il + 2Kx)GDiinJ) + ^GD2in,l) 



■m. 



(A.ll) 



where Guiiji^l) and GD2{n,l) are given as, 

r ( n-f'^V' ^ -^W2{n-i-j-iy 

\ 21+3 



r(4/ + 2 J + 3) 



\narJ n22(^ 



Kp \ m 



l-kj 2 



r(2/ + 3)| 



nttr 



2 + mnar 



F -2j, 21 + 3; 4/ + 3; 



2 + mnar 



m^k'^tvn 



nar \ 



21+2 



-T{21 + 2) ( _- ) F ( -2j, 2/ + 2; 4/ + 3; 
^1 — k'^Y \2 + mnar/ V 2 + mnar 



m?k'^Kp 
il-k^V 



r(2/ + 2) 



nar 



2 + mknar 



21+2 / 4 

F -2j, 21 + 2; 41 + 3; 



2 + mknar 



GD2{n,l) 



f 2 



V rictr / n 



2' 1 "^1 f2{n - / - J - 1)\ r(4/ + 2j + 3) 

22("-'-i) ^t^o 1^ n_/_j_i jj!r(2/+j + 2)r(4/ + 3) 



X 



na^ 






2Z+3 



F -2j, 2/ + 3; 4/ + 3; 



2 + mna^ 
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m^k^htp „,,, ^, / nttr \2«+2 



T(2l + 2) ( ^^^ I F 1—2? 2/ + 2" 4Z + 3' 
(1 — A;^)^ v2 + ?TTOar/ \ ' ' ' 2 + rrina^ 

I ^'^'^P r(2f + 2) f ^"^ ) ''^' F ("-2^, 2/ + 2; 4/ + 3; ^ ) 

(1 — fc^)^ \2 + mknarJ \ ' ' ^2 + mknar/ 

Similarly, the fine structure potential with form factors gives rise to flTSj) where, 

^ 2 ^^^+3 1 "^V2(n-/-j-l)\ r(4/ + 2j + 3) 



Knar J 2ri22("~'-i) ^^^ \^ n-l-j-1 J j\T{2l +j+ 2)1(4/ + 3) 



1 — k'^Y j \2 + mnar J V ' ' ' 2 + mna 

-r(2/ + l)mfl + -^)(-^^^^)''^'F(-2j,2/ + l;4/ + 3;— i 

y (1 — K^j^y V2 + mna.r/ V 2 + innar 

.m^ f K„ \ / nar \^'+^ / , , 4 

-r(2/ + 2)— 1 + -^ F -2j, 21 + 2; 4/ + 3; 



2 V 1 — k"^ J \2 + mnar J V ' ' ' 2 + mnar 

+r(207^^bT^h^— f— F -2j,2/;4/ + 3; 



1 — k'^y \2 + mknar J \ ' ' ' 2 + rri/cna^ 



r\iT) I I l/iAiFf 



+T{2l + l)mk- — ^ r ^ F -2j,2/ + l;4/ + 3 



[1 — k'^Y \2 + mknar J \ ' ' ' 2 + mkna 



. (A.12) 
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